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ABSTRACT 

We examine the improvements in constraints on the Hnear growth factor G and its 
derivative / = dhiG/dlna that are available from the combination of a large-scale 
galaxy redshift survey with a weak gravitational Icnsing survey of background sources. 
In the linear perturbation theory limit, the bias- modulation method of McDonald & 
Seljak allows one to distinguish the real-space galaxy clustering from the peculiar 
velocity signal in each Fourier mode. The ratio of lensing signal to galaxy clustering 
in transverse modes yields the bias factor b of each galaxy subset (as per Pen 2004), 
hence calibrating the conversion from galaxy real-space density to matter density in 
every mode. In combination these techniques permit measure of the growth rate / 
in each Fourier mode. In principle this yields a measure of the growth rate that is 
free of sample variance, i.e. the uncertainty in / can be reduced without bound by 
increasing the number of redshifts obtained within a finite survey volume. In practice, 
the gain from the absence of sample variance is bounded by the limited range of bias 
modulation among dark-matter halos. Nonetheless, the addition of background weak 
lensing data to a redshift survey increases information on G and / by an amount 
equivalent to a 10-fold increase in the volume of a standard redshift-space distortion 
measurement — if the lensing signal can be measured to sub-percent accuracy. This 
argues that a combined lensing and redshift survey over a common low-redshift volume 
of the Universe is a more powerful test of general relativity than an isolated redshift 
survey over larger volume at high redshift, especially as surveys begin to cover most 
of the available sky. An example case is that a survey of « 10^ halo redshifts over half 
the sky in the redshift range z = 0.5 ± 0.05 can determine the growth exponent 7 for 
the model / = il^ to an accuracy of ±0.015, using only modes with k < 0.1ft- Mpc^^, 
but only if a weak lensing survey is conducted in concert. 

Key words: gravitational lensing, methods: statistical, large-scale structure of Uni- 
verse 



1 INTRODUCTION 

The growth of large-scale structure in the Universe is a 
competition between gravitational attraction and the expan- 
sion of the Universe. In the linear-perturbation-theory limit 
(LPT), each Fourier mode evolves independently, and the 
mass density fluctuation amplitude at a comoving wavenum- 
ber k is (5(k, a) — G{a)5{'k,ao), with a the scale factor. 
Knowledge of the linear growth function G{a) would, in 
standard general relativity, reveal the expansion history of 
the Universe, including any efl'ects of "dark energy." Given 
prior knowledge of the expansion history a{t), the measure 
of G{a) tests the time and scale dependence of any devi- 
ation from General RelativityQ The surprising observation 



that the growth a{t) is accelerating in recent epoch strongly 
motivates tests of the accuracy of General Relativity at the 
largest observable physical scales, preferably to accuracy 
< 1%. 

Since the growth of large-scale structures is unmeasur- 
ably small over the brief history of human astronomical ob- 
servations, we cannot directly track the growth in a given 
mode. Growth measures must adopt one of two alternative 
strategies. The first is to compare fluctuations at the same 
k at two different epochs on the light cone, invoking the 
Cosmological Principle to permit comparison of structure 
at distinct locations. The comparison is therefore limited by 
the sample variance from the finite number N,n of mode am- 
plitudes that are drawn from the distribution and observ- 
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^ We will ignore possible scale dependence of G{a, k), that could 
arise from modifications to General Relativity, but this is merely a 



notational convenience: all the results in this paper apply equally 
well to a scale-dependent G or /. 



2 Bernstein & Cai 



able at each epoch: the power P is measured to accuracy 
ap/P ^ y^2/Nm regardless of the noise level in the ob- 
served fluctuation field. If the theoretical prediction for P is 
sufficiently accurate for k < fcmax, then the number of modes 
available to a survey of volume V is Nm = V^fcmax/Svr'^. 

The sample variance from the finite observable volume 
of the Universe is often called cosmic variance and places 
a fundamental limit to the measure of growth by power- 
spectrum comparison in a chosen k range. The sample vari- 
ance limits become particularly acute at low redshifts, be- 
cause the available volume scales as z^Az and because non- 
linear growth of structure reduces the fcmax at which accu- 
rate theory is possible. This is unfortunate since the mani- 
festations of dark energy or modified gravity are thought to 
be the largest at the present epoch. 

The second strategy for measuring growth is to invoke 
the continuity equation to link the growth rate to the veloc- 
ity amplitude in each mode. If each mode evolves indepen- 
dently, and the mass is conserved with a single-valued pecu- 
liar velocity v at each location, then the continuity equation 
^-|-V - (l-|-5)v = becomes, to first order in the density 
perturbation, V • v = — ^ . In Fourier space, it yields 

v(k) = ///5(k)g, (1) 

where / = g'"^ is the growth rate and H is the Hubble 
parameter at that time. In real space or Fourier space, the 
velocity field is proportional to the mass density field, with 
the constant of proportionality / delivering the desired in- 
formation on the (differential) growth of structure. Since 
the LPT velocity field is irrotational, it suffices to map one 
component of the velocity field. Massive observational ef- 
forts were undertaken in the 1980's and 1990's to map the 
density and (radial) velocity fields of galaxies in the nearby 
(< lOO/i"^ Mpc) Universe. These were stymied by (among 
other issues) two difficulties: first, the observed radial ve- 
locity of a galaxy at distance r is Voba = v ■ r/\r\ + Hr, so 
these surveys required distance indicators for each galaxy to 
determine r and the (radial) peculiar velocity v indendently, 
introducing large statistical and systematic errors. Second, 
even with perfect distance indicators, the density field would 
be derived from the galaxy field, so the ratio of velocity to 
density will yield f /h, not /. Inferences on the growth rate 
/ can be only as accurate as prior knowledge of the galaxy 
bias factor b. 

Current proposals for measuring / avoid the need for 
distance indicators — thus extending the measurements to 
larger volumes and higher redshifts — by analyzing the direc- 
tional dependence of the redshift-space galaxy power spec- 
trum P''(k). If we survey a tr acer populati on that follows 
the velocity field of the matter. [Kaiseij (| 19871 ) shows that the 
LPT value of its Fourier fiuctuation amplitude tJ" in redshift 
space is related to the real-space amplitude 5^ of the tracers 
and 5 of the mass by 

5^(k) = <5'-(k) + /M^5(k). (2) 

Here ^ is the cosine of the angle between k and the line of 
sight (taking the plane-parallel approximation). The stan- 
dard simple approach is to next assume that the tracer den- 
sity has some bias b with respect to the matter, so that 
5^ = b5, but faithfully traces the matter velocity field on 
large scales, in which case the redshift-space power spec- 



trum becomes 

p=(k) = {b' + 2bf^J? + fi,^)P{k) 

= fp{k)[{b/ff + 2{b/f)^^ + ^,^]. (3) 

Measurement of Ps{k,fj,) again cannot constrain / without 
prior knowledge of b, however the combination f^P can be 
constrained. Current state of the art is demonstrated by 
the ~ 20% measures of f^P in bins of width Az — 0.2 by 
(|Blake fc et al.ll201ll ). We have fP = fG^Po, where Po is 
an "initial" power spectrum at an early epoch, e.g. recombi- 
nation, and G is the linear growth since that epoch. Relying 
on the CMB to determine Po, we can constrain the growth 
quantity fG by comparing the power in distinct modes (dif- 
ferent fj, instead of different epochs this time), which again 
leads to a fundamental sample- variance fioor. Note also that 
any uncertainty in Po, for example from uncertainty in the 
reionization history, propagates into an error of half the size 
in fG. The regression against fi necessary to extract fG 
amplifies the sample variance significantly. A Fisher anal- 
ysis of the standard RSD ana lysis in the sample-variance 
limit yields ain/G ~ ^21/Nm for the case b/f — 1.4 (with 
the prefactor becoming even le ss favorable for larger b/f). 

iMcDonald fc SeliakI (|2009l )[MS] propose an important 
improvement to the peculiar- velocity measurement. Divide 
the galaxy sample into subsets with differing bias factors bi, 
and assume that the overdensity of each bin in redshift space 
follows the Kaiser form: 

St = (b, + ftiy + e,, (4) 

with (eitj) = /rti describing a diagonal stochasticity ma- 
trix. For the usual assumption of Poisson sampling, rii is the 
space density of the objects in bin i. In the limit of l/n; — >■ 0, 
MS note that observation of the 51 in a single mode ai ji = Q 
will yield the values of bi up to an unknown normalization 
b. Then a subsequent noiseless observation of another single 
mode at /i 7^ will allow one to regress 5| against b to de- 
termine both f^J?5 and 65, as illustrated in Figure [1] Since 
is known, in the limit of low stochasticity this method 
produces 

(i) a measure of f^P — f^G^Po yielding uncertainty 
'^in fG = l/V^Nm, assuming known Po. This is the sample- 
variance limit if we could simply view the "naked" velocity 
field distinct from the real-space clustering, and is a factor 
of ~ 40 improvement in variance over the standard (RSD) 
analysis, which was degraded by the need to regress P" 
against /i; 

(ii) a measure of f /b from taking the ratio of velocity 
to real-space density amplitudes in each mode. Since the 
measure can be made mode-by-mode, it is not limited by 
sample variance, i.e. can in theory reach unlimited precision 
from a finite set of modes. 

When the Ui are finite, these gains are of course ameliorated. 
Numerical results are published for constraints on fG ex- 
pected from values of Ui and bi that are ar bitrarily assigned 
in 2 b ins (MS) ; or for halos split into 2 bins (|Gil-Marm et al.l 
I2OIOI): or for pote ntial observed samples split into 4 bins 
l|White et al.ll2009l l. In this paper we derive the optimal re- 
sults attainable for the MS method using any number of bins 
of halos. 

The sample- variance-free measure of f/b is not useful 
for cosmological constraint without further information on b. 
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Figure 1. Schematic illustration of the McDonald-Seljak tech- 
nique: in a chosen Fourier mode with real-space mass density fluc- 
tuation (5, we observe the redshift-space amplitude 5^ of galaxies 
with different biases 6;. If the Kaiser formula ^ holds, then lin- 
ear regression of the S? data points against bias will yield the 
y-intercept value flJ,^5 and the x-intercept value //i^. The former 
is one sample from a Gaussian with variance fi^P. The latter 
gives / without sample variance. However the uncertainty in / 
from this mode's data is amplified if ^ ^ 1, or if the range Afe of 
galaxy biases is small compared to the typical b + f value. 



Cross-correlation of galaxy redshift surveys with weak grav- 
itational lensing shear fields — which are a direct measure of 
the true mass overdensity 5 — yields a sample-variance-free 
measure of 6 (|Penll2004l ). and MS comment that combining 
this technique with the multiple-tracer RSD measure could 
yield improved cosmological constraints. The main result 
of this paper is to quantify the gain in growth-of-structure 
constraints from this combination of galaxy redshift surveys 
with weak lensing surveys. This combination is attractive 
for several reasons: 

(i) The combination of bias-modulated RSD data with 
lensing data yields an estimate of / that is free of sample 
variance, meaning arbitrarily good accuracy with finite sur- 
vey volume. This is very attractive for measuring growth on 
very large scales at low z, where Fourier modes are scarce. 

(ii) This combination measures / directly, breaking de- 
generacies with b or P. Consequently the cosmological model 
can be constrained without being subject to uncertainties in 
the CMB power spectrum normalization due to reionization 
complexities, potentially even exceeding the cosmic variance 
level of the CMB measurement. 

(iii) Breaking the degeneracy between / and P enables 
a consistency check between the growth function G and its 
derivative / over the observed redshift range. This allows 
some of the theoretical assumptions to be checked, for ex- 
ample the assumption that the continuity equation for mass 
can be used to infer the velocity field for halos or galaxies. 

We first give a qualitative discussion of the types of mea- 



surements that are free of sample variance, then we provide 
the analytic formalism and derive the scaling of the errors 
on / and P in the case of independently-determined biases. 
We do not limit our analysis to 2 bins of galaxy bias. We 
next provide a numerical analysis of the growth constraints 
available in a candidate redshift survey of all halos above 
mass Mmin in a volume centered at z = 0.5. This is a red- 
shift at which the standard RSD analysis is becoming very 
constrained by cosmic variance. Unless noted otherwise, we 
assume a fiducial flat ACDM cosmology with the following 
parameters: n„=0.27, = 0.73, fib = 0.046, as = 0.9, 
ns = 1, Ho^ 72. 



2 SAMPLE-VARIANCE-FREE 
MEASUREMENTS 

Consider a random variable S that is drawn from a zero- 
mean Gaussian distribution with variance P. If we draw 
{5i,52, ■ ■ ■ ,5ns} from the distribution, and execute a mea- 
surement process on each sample that yields Si — 5i + due 
to some measurement noise e^, then how well can we deter- 
mine the power P? If the noise power £ = (e^) is known, 
then the uncertainty on P satisfies 



2 



{i + £/py 



(5) 



The constant term inside the parentheses sets a sample- 
variance limit that cannot be reduced by better measure- 
ments, only by obtaining more samples. If 5 is the mass den- 
sity or gravitationl potential fluctuation in a Fourier mode, 
then Ns is the number of independent modes in the surveyed 
volume, which we call Nm'- 



V A-Kki, 



(27r) 



(6) 



Next consider the case where two tracers of the under- 
lying random variable 5 are available. Two measurements, 
— Si + ef and 5^ = bSi + ef , are available for sample 
i. Tracer B is biased by some factor b with respect to the 
underlying field S, and measurement errors are uncorrelated. 
The covariance matrix between 5"^ and 5^ is 
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P + £^ bP 
bP b'^P + S'^ 



(7) 



If we assume that is known but that P, 6, and £^ are 
free parameters, we can derive the covariance matrix for our 
estimates of these parameters using the Fisher matrix for- 
malism. The covariance matrix for these three parameters 
is, at bes t, the inverse of the Fisher matrix F having com- 
ponents (|Tegmark et al.lll997^ 



— -^^"C iC A^C AB aC a\C A 



B.J 



(8) 



The uncertainty in InP in this two-tracer case remains 
the same as (O, and has the same sample- variance limit 
finP ^ \/ 2 /Ns ■ The uncertainty on the bias parameter b 
(after marginalization over P and is 



2 



1 



cA cB 

-p+l^ + ois/py 



(9) 
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Note that the leading order term in £/P is Unear: the 
sample- variance term is missing, and a decrease in measure- 
ment error reduces uncertainty in b without bound. The fun- 
damental difference between b and P is that the former can 
be measured accurately from a single pair of samples 6^ and 
5^ whereas the latter is a population property. An intuitive 
result is that, with no measurement noise, b = 5^/5^ can be 
measured perfectly as the slope of the locus of measurements 
in the — 6^ plane. 

We take advantage in this paper of three circumstances 
in which cosmological information can be extracted from 
the relative signals of distinct tracers of the mass density 
fluctuation field. 



2.1 Weak lensing convergence vs galaxy 
distribution 

For fluctuation modes transverse to the line of sight, gravita- 
tional lensing imparts a shear fleld on background galaxies 
that can be measured with weak lensing techniques. The 
associated lensing convergence amplitude k is equal to the 
matter fluctuation 5, modulo some distance factors, which 
we will assume are known. The second mass tracer is the 
galaxy density map, wh ich has an unkn own bias factor b 
with respect to matter. IPenI (|2004h and iBernstein fc JainI 
l|2004 ) both note that the proportionality constant between 
lensing and galaxy signals can be measured without sam- 
ple variance. For the lensing signal, the measurement er- 
ror £^ is from shape noise, and at large scales {£ « 100) 
the ra tio £/P for "cosmic shear" power is < 0.1 in each 
mode jTakada fc Jainll2Q04l '). For the galaxy tracer, the mea- 
surement error £^ arises from stochasticity in the galaxy 
distribution, often assumed to be Poisson noise, in which 
case £^/b^P = (n&2p)-\ where n is the space den sity of 
galaxies, and this can be quite small. ICai et all (|201ll ) show 
that an optimally-weighted survey of halos can reproduce 
the mass density with very low stochasticity — e.g. an opti- 
mally weighted survey of halos with M > 10^'''^h~^ A'Iq can 
obtain £/b^P < 0.1. Equation (O then suggests that the 
comparison of lensing and galaxy information in the ~ 10* 
transverse modes with I ^ 100 should yield bias uncertain- 
ties (Tlni, < 0.01. 

We defer to future work the examination of the com- 
plexities in constraint of galaxy bias by this method, such 
as the simultaneous solution for bias and distance factors, 
and the superposition of many galaxy redshifts along a given 
lensing line of sight. In this work, we will simply assume 
that large-scale galaxy biases — or at least the bias of a cho- 
sen weighted combination of halos — can be derived by the 
lensing+galaxy surveys to potentially sub-percent accuracy, 
limited by signal-to-noise rather than by sample variance. 

2.2 Real-space clustering for different bias bins 

Consider next a redshift survey that has targets divided into 
bins expected have different bias, as per MS. We will con- 
sider, optimistically, that redshift targets are dark-matter 
halos of known mass, and that we bin the targets by halo 
mass. For a transverse mode = 0), the redshift-space sur- 
vey has no velocity contribution [Equation Q]. As noted 
by MS, the ratio bj/bi of biases of two different mass bins i 



and j can be measured without sample variance, since this 
is the ratio of two tracers of the same underlying field 5. 
Hence the biases of all halo bins can be determined, without 
sample variance, up to an overall scaling factor b. The ad- 
dition of lensing data allows a high-precision determination 
of b as described above, again without sample variance, so 
the combination of both techniques will yield biases of all 
halo-mass bins. 



2.3 Velocity field vs density field 

MS offer a third means of extracting cosmological informa- 
tion from a ratio of tracers. With the redshift-survey target 
galaxies divided into bins of differing bias, it is possible to 
distinguish the two terms bi5 and f^J.^5 that contribute to 
the amplitude of the redshift-space amplitude 5f in every 
mode. The ratio of these two components is ffi^/bi, another 
quantity that can therefore be measured without a sample- 
variance limitation. If the bi are known, / is obtained. 

We see therefore, that combining a lensing survey with 
a redshift survey over the same volume provides a means to 
measure the ratio of the velocity signal to the mass over- 
density in every (non-transverse) mode. This ratio gives /, 
without any fundamental limit from sample variance, with- 
out the need to compare to a power spectrum at a reference 
epoch, and without any dependence on bias. 



3 FORECAST METHODS 

We construct a Fisher information matrix for the constraint 
of growth-related parameters by a canonical redshift survey 
of Fourier modes with k < fcmax in volume V at redshift z. 
The redshift-space structure is taken to follow the Kaiser for- 
mula (U), and the real-space matter fluctuations are taken to 
obey {(5^(k)) = G^Po(k) with no correlation between modes. 
The shape Po (k) of the power spectrum is taken as known, 
while the amplitude G, i.e. the growth function, is a pa- 
rameter of interest. The fluctuations in a single Fourier 
mode of the survey volume are distributed, in LPT, by a 
multivariate zero-mean Gaussian with covariance matrix 

= Cov{St,S'^) = (6. + f^^)(bj + fy?)G'Po + £,,(10) 
£ij = {€iej) (11) 

The model has free parameters (which we index by Greek 
letters) 9a £ {/, G, 6i, 62, • ■ • , &]v}. In the model where the 
spectroscopic targets are halos that sample the mass field by 
a pure Poisson process, we will have the stochasticity matrix 
£ — diag(fi) = diag(l/wi), where rij is th e space density o f 
targets in bin i. iHamaus et al.l (|2010h and ICai et al.1 (|201lh 
demonstrate that the simple diagonal Poisson-noise formula 
for £ does not fully describe the stochasticity of halos in A'^- 
body simulations. We have verified that one can treat the 
diagonal elements £i as free parameters of the model, with 
insignificant degradation of the growth measurement, so this 
part of the Poisson assumption is not critical. The success 
of the analysis described herein does, however, depend cru- 
cially on the assumption that the off-diagonal elements of £ 
vanish, so future work will need to quantify the corrections 
needed for non-diagonal stochasticity. 

The Fisher matrix for information from a single mode 
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follows th e standard form for m ultivariate, zero-mean Gaus- 
sian data jTegmark et al.|[l997l '): 



a0 



-Tr 



i9C i9C 



(12) 



The total Fisher matrix is then the sum over all independent 
k modes in the survey volume. 

The total Fisher matrix is then used to predict the co- 
variance matrix of In / and In G after marginalization over 
all other parameters, namely the biases. Marginalization and 
mapping to new parameters are done with s t andar d tech- 
niques, e.g. as summarized bv lAlbrecht et all (|2009l) . 



3.1 Cases 

We will forecast uncertainties on the growth parameters for 
the following types of experiments: 

• Standard RSD: All halos with mass M > Afmin in 
the survey volume have redshift measurements and are com- 
bined into a single sample. We marginalize over the bias of 
this combined sample. We know that there will be a degen- 
eracy between In / and In G since only the quantity f^P can 
be constrained in this experiment. We will report the uncer- 
tainty on the quantity In /G. As noted in the int roduct ion, 
there is a sample variance limitation crin/c ^ y^2l/Ns for 
b/f = 1.4. 

• MS: The same as Standard RSD, except that halos 
are binned by mass and the analysis suggested by MS is 
performed. Biases of every mass bin are free parameters. 
We typically use 10-20 mass bins, which will produce results 
indistinguishable from the ideal limit of an analysis without 
information loss due to binning. Again In/ and InG are 
degenerate without prior knowledge of biases, so we report 
uncertainties on In/G. The sample variance limitation is 

• Fixed Biases: The same as the MS analysis except this 
time we assume that biases of all halo mass bins are known 
perfectly through combination with a weak lensing survey 
of the same sky area. In this case we can simply strike the 
bicis-related elements from the Fisher matrix. The degener- 
acy between / and G is broken so we can report two un- 
certainties (Tin/ and (TinQ. Constraint on alternative growth 
scenarios must account for possible covariance between In/ 
and InG in the measurement as well. The sample- variance 
limit for ainc is 1/ V^Nm, the same as for /G in the MS 
method. 

• Weighted Bias: A more realistic scenario than the 
Fixed Bias case is that the weak lensing measurements pro- 
vide an estimate of the bias of a single weighted combina- 
tion of the halos. This would be found by cross-correlation 
of the lensing convergence with a weighted map of halos. 
We assume that the bias of this weighted halo population 
is determined to a fractional accuracy of aint- The effect on 
the Fisher matrix is to add WiWjbibj/af^f, to the element 
of the Fisher matrix linking bias bi with bias bj. Here Wi is 
the weight on halos in bin i, normalized so that Y^Wi = 1.. 
We marginalize over all the biases after adding this prior 
on the weighted bias. We assume that mass weighting 
would be used, as a simple appro ximation to th e minim um- 
stochasticity weights derived by iHamaus etall (|2010l ) and 



ICai et al.l (|20lj ). The weighted-bias-measurement experi- 
ment breaks the degeneracy between / and G. 

3.2 Algebraic Results 

Algebraic expressions for the Fisher and covariance matrices 
in these cases are available, because the matrix G is invert- 
ible via the Sherman-Morrison formula. Defining a vector u 
with Ui — bi + ffi^, we have 

T 



C 



= £ 



£ + uG^Pou 

G^Po£-^uu^£-^ 



(13) 
(14) 



l + G^PouT£'^u 

In the Poisson case, £~^ — diag(?ii) and the Fisher-matrix 
elements all can be expressed in terms of number-weighted 
sums over the halo population. These results are stable un- 
der change in bin size as long as the bias values do not vary 
widely within any given bin. Marginalizing over bias and 
summing over the distribution of /i for all the modes leads 
to rather lengthy, opaque expressions for the final uncertain- 
ties. However in the case of a single mode at fixed /i with 
known biases, the uncertainties in G and / can be derived 
in an illuminating form. If we define n = ^ rii as the total 
space density of halos surveyed, we obtain 



2 



l + nP{{b + f^l^f) 
^i^{nPYY&r{b) 

1 {(b+ff^r) 



(15) 



Var(6) 
(nP> 1). 



(nP> 1), 



(16) 
(17) 



fi^nP 
1 

G2 2 

The averages and variance in this expression are taken over 
the targets of the redshift survey. 

This scaling is understood from examination of Fig- 
ure [Jl illustrating the MS method. The quantity of interest, 
//i^, is equal to the 6-intercept of the linear regression of 
the fiuctuation amplitudes Si against the bias values hi. The 
main uncertainty A{ffi^) in the location of this intercept 
come from statistical fluctuations Am in the fitted slope m 
of 5i vs bi: 



A(/m') = t^'Af . 



Am. 



(18) 



The uncertainty Am in the slope will be roughly (noise in 
5i)/(span of observed b values), or (Am)^ ~ [nVar(6)]~^, 
since the Poisson variance in Si is 1/n. Also we note that 
(m^) = P since S is the slope of the line. Putting these 
together. 



M^(A/)= 



1 {b+ff.y 



(19) 



nP Var(b) ' 

which is very similar to 1)17^ . We see that a small value of 
Var(fe) in the target population will lead to a large lever-arm 
in determining the 6- intercept that defines /, hence a narrow 
range of biases will produce poor constraints on /. 

Equation (|17p also exhibits the expected characteris- 
tic of a sample- variance- free measurement, namely that af 
drops without bound as the measurement noise is driven to 
zero (n — oo). We will see, however, that this gain is not 
realized, because a real halo population has Var(b) — >■ as 
we seek large n by going to ever-lower halo masses. 
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The same single-mode analysis shows that the limit 
nP — >■ oo yields an uncertainty in the power spectrum am- 
plitude ainG — >■ l/%/2, i.e. a sample variance limit equal to 
that for In fG in the MS method. 



4 RESULTS 

Figures[2]plots the uncertainty in the growth measures vs the 
minimum mass Mmin of halo included in a fiducial redshift 
survey. We make the following assumptions for the fiducial 
survey: 

• The survey covers /sky ~ 0.5 of the celestial sphere and 
a redshift range Az — 0.1 centered at z = 0.5. In the fiducial 
ACDM cosmology, the survey volume is 2.6 h~'^Gpc^, and 
/ = 0.72 at the survey redshift. 

• The survey obtains a center-of-mass redshift for all 
halos with M > Mmin. The number and bias of ha- 
los v s mass are t aken from a Sheth-Tormen mass func- 
tion jSheth fc Torm cn 1999;), with model parameters listed 
therein. 

• The matter fluctuations have a power spectrum given 
by linear perturbation theor y, with the transfer functi on 
computed from CMBFAST (|Seliak fc Zaldarriaga|[l99^ V 

• The halos populate the matter distribution by a biased 
Poisson process. 

• The halos follow the LPT matter velocity field without 
additional bias or stochasticity, so the Kaiser formulae hold. 

• Only modes with k < 0.1/iMpc~^ (optimistic) or A; < 
0.03/iMpc~^ (pessimistic) have theoretical redshift-space 
predictions of sufBcient accuracy to be useful for percent- 
level constraints of growth. A^-body simulations show 20% 
departures from the Kaiser formula for redshift-space distor - 
tions of halos at < 0.1/iMpc"^ (|Okumura fc Jind [20T1I ) . 
with significant dependence on halo mass. Substantial ad- 
vances in theoretical understanding of nonlinear redshift- 
space effects will be needed to extract percent-level growth 
constraints at these scales fe.g. Ijennings et al.ll201ll ). 

It is desirable to measure the growth function to percent- 
level accuracy in several redshift bins of this width at 2 < 
1 in order to test general relativity (in combination with 
percent-level constraints on a{t) from distance-measurement 
methods). We note that the total number of independent 
Fourier modes in the fiducial survey volume is only Nm = 
1150 (43,000) for fc^ax = O.OS/iMpc"^ (0.1). We recaU that 
the standar d RSD method has a sample- variance floor of 
o"in/G ~ y^21/Nm and therefore a 1% measure of fG is 
not attainable for fcmax ^ 0.1/iMpc^^. The sam ple variance 
limit of the MS method is a\nfG = \/ 1 /'^Nm and the 1% 
meaurement is potentially achievable in the optimistic case. 

4.1 Surveys without lensing 

Figures [2] plot the uncertainty in In fG for the Standard 
RSD and MS methods, and then plot ainf and aino inde- 
pendently for the Fixed Bias and the Single Bias Measure- 
ment cases. In this last case, we plot results from two levels 
of uncertainty on the weighted bias b. Two horizontal axes 
plot the Mmin for the halo survey, and the total number of 
halos in the volume, i.e. the number of redshift measure- 
ments required for the survey. 



The Standard RSD analysis (in blue) approaches a 
sample-variance-limited plateau as expected, with very lit- 
tle improvement once Nh > 10^-10*^ halo redshifts are ob- 
tained. 

The MS analysis (purple) shows significant improve- 
ment over the Standard RSD analysis when > 10^ redshifts 
are obtained, or halos < Mq. The cosmic variance 

limit for the MS analysis is fin/c > l/V^Nm (dotted pur- 
ple line) but this is not attained even with surveys of > 10® 
redshifts. This is attributable to the limited number of halos 
with 6 > 1, as explained below. In the conservative case of 
fcmax = 0.03/1 Mpc~^, uncertainties of < 10% in fG require 
Nh > 10^ redshifts to be measured in this bin for the MS 
method. If the theory supports fcmax ~ O.l/iMpc" , then 
only Nh < lO'' redshifts are needed to reach 10% accuracy, 

and 10^ redshifts yield a 2% measure of fG. 

1/2 

Uncertainties on fG scale roughly as A'^^ when using 
the standard or MS methods at A^^ < 10^. In practice a 
survey has a trade between depth and sky coverage, and 
one can ask whether better cosmological constraints result 
from a large, sparse redshift survey or a deeper, smaller-area 
survey given a fixed survey cost (duration) . If every redshift 
can be attained at equal cost, then this scaling implies a 
flat trade, i.e. constraint on fG are independent of survey 
area when Nh < 10^. However the a-vs-Nh curve flattens 
at higher Nh, which indicates that one should favor area 
over depth when Nh > 10^ . This would be even more true if 
redshifts are more expensive to obtain in lower-mass halos. 



4.2 Surveys with lensing data 

The red lines show the uncertainties on In G (solid) and In / 
(dashed) in the Fixed Bias case of perfect lensing calibra- 
tion of all galaxy biases. We note first that constraints on 
In / (dashed red) are obtained at accuracy similar to the 
In fG constraints with the MS method at equal Nh . Next 
we note that the In G constraint (solid red) approaches a 
sample-variance limited plateau much more rapidly than 
does the MS method. For fcmax = 0.03/iMpc^^, one can ob- 
tain 3% errors in In G with a survey of only ~ 10'' halos at 
M > 10^*/i~^Mq, 10* times fewer redshifts than are needed 
for similar fG constraints using the MS method alone, and 
completely impossible with Standard RSD methods. This is 
attributable to the lensing information breaking the degen- 
eracy between bias and matter power spectrum amplitude. 

The dashed red line shows that the constraint on / does 
indeed continue to decrease without apparent bound as Nh 
increases and shot noise is reduced. It is disappointing, how- 
ever, that (Tin/ is larger than ainc until the survey exceeds 
10^ redshifts in this bin. Using Equation p7p as a guide to 
the behavior of ain / , we can see that the survey can attain 
very large nP but still have the / constrain t degra ded by 
an unfavorable lever arm factor //i2)/^Var(6). While 
massive halos have 6 > 1 and a significant range of bias, all 
of the halos with masses IOIO-IO12/1-IM0 have b within a 
few percent of unity. As low-mass halos dominate the full 
sample, Var(&) drops, degrading the gains from decreased 
shot noise. The scaling is roughly ainf oc N^^^'"^. 

The green lines show the effect of degrading perfect 
knowledge of all biases to finite errors on a single weighted 
bias combination. We find that the measure of a single 
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Figure 2. Forecasted uncertainty in growtli parameters are plotted against the number of halo redshifts obtained in our fiducial survey: 
z = 0.5, Az = 0.1, /sky = 0.5. The survey is assumed to reach all halos above a minimum mass, marked on the top axis. Left and 
right panels assume difl'erent wavenumbers fcmax to which linear perturbation theory is sufficiently accurate. The legend marks the 
different types of analyses used to extract redshift-space distortion information: blue and purple curves use single-bin or McDonald- 
Seljak bias binning to analyze pure galaxy-redshift information and extract the degenerate combination fG. With weak-lensing enables 
measurements of galaxy bias, / and G can be measured separately to the plotted accuracy, for perfectly known biases (red) and for a 
single weighted bias uncertain to the marked levels (green). Combination of weak lensing data with the galaxy redshift survey not only 
enables a direct measure of /, but a distinct measure of G with substantially better precision. The upper horizontal dotted lines mark the 
cosmic- variance limits for the standard RSD measure, and the lower dotted lines are the sample variance limits for In fG (MS method) 
or In G (lensing-l-redshift methods) . 



weighted galaxy bias results in growtli constraints that are 
nearly as good as the Fixed Bias case. The degradation in 
accuracy of parameter 6 G {/, G} is insignificant as long as 
o"6 -C crin9. In particular, this means that the weighted bias 
should be measured to 0.01 (0.001) accuracy for fcmax = 0.03 
(0.1). This is not a trivial requirement — a future paper will 
detail the attainable accuracies in b from the lensing exper- 
iments. 



4.3 Modified gravity constraints 

How do we compare the merits of constraints of / to con- 
straints on G or on /G? For the lensing+galaxy surveys, 
does the correlation between measures of / and G enhance 
or degrade the detection of modifcations to General Rela- 
tivity (GR)? We address both questions by quantifying the 
constraints on a simplistic model fo r deviations from GR . 
iPeeblej (|l980l ): iLahav et all ([l993); iLinder fc Caiml (|2007l . 
and references therein) note that the growth equation for 
GR is very well approximated by 

f = ^l. (20) 

with 7 in a narrow range near 0.55 for a wide variety of 
dark-energy models. Gravity theories other than GR tend 



to have a very different value of 7. Therefore, a constraint 
on 7 may provides sensitive test of GR. 

A different choice of 7 with fixed ACDM expansion his- 
tory would alter both G and / in a predictable way. We find 
numerically that d{ln f) / d'y ~ 3d{\nG)/d'y in this model, 
so we can genetically expect constraints on the differential 
growth / to be 3x more valuable for detecting deviations 
from GR than constraints on the integrated growth G and 
3/4 as valuable as fG constraints of equal fractional ac- 
curacy. Figures [3] plot the uncertainty on 7 that would be 
implied by the forecasted / and G measurements from our 
fiducial surveys. 

The plots suggest that a coincident weak lensing survey 
greatly enhances the ability of redshift surveys to constrain 
deviations from GR: the lensing determination of a weighted 
galaxy bias to 0.01 (0.001) accuracy yields a cr^ value that is 
at least vTOx better than the standard RSD constraints for 
a survey of fixed fcmax = 0.03/i Mpc^^ (0.1) and fixed num- 
ber of redshifts. In other words, adding a coincident lensing 
survey to a standard RSD survey improves its growth con- 
straints by an amount equivalent to a 10 x expansion of the 
sky coverage of the RSD survey. A redshift-only survey us- 
ing the MS technique can recover some of the advantage of 
the combined survey, but rather slowly: the MS approach 
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Figure 3. Forecasted constraints on deviations of the growth parameter 7 in Equation l|20|l from the General Relativity value of 7 = 0.55. 
The axes and line colors apply to the same fiducial survey and types of analyses as in Figure [2j In the case of fcmax = 0.03hMpc-i (0.1), 
a lensing-based bias measurement to accuracy 0.01 (0.001) provides improves constraints on 7 over a standard redshift-space-distortion 
experiment by an amount equivalent to a 10-fold increase in survey size. 



requires 10*-10^ redshift measurements to attain the 7 pre- 
cision that a combined survey reaches with 10® redshifts. 

The scahng of with survey depth is instructive. For 
shallow surveys (M^in > lO"/i"^M0, iV^ = lO^-lO'^), the 

— 1/2 

combined survey gains precision quickly, roughly as A^^ , 
the slope at which the depth-vs-area tradeoff will be rather 
flat. The gains in precision for deeper combined surveys are 
quite shallow, aj oc Nj^^''^ . This reduced return on invest- 
ment is traceable to the small variation in bias for lower- 
mass halos. Note that the inflection occurs for a survey that 
requires redshifts only for galaxy-cluster-mass halos, which 
are easily identified with imaging surveys at z < 1. If biases 
of such halos can be determined to 0.01 (0.001) fractional 
accuracy from a coincident lensing survey, then uncertain- 
ties in 7 of ~ 0.1 (0.015) per Az = 0.1 bin can be obtained 
with only ~ 10^ redshift measurements per bin. Both results 
would require > 10* redshifts per bin if using the MS tech- 
nique without lensing, and are better than would be possible 
from any standard RSD survey, regardless of depth. 



5 DISCUSSION 

The combination of bias-modulated redshift-space measure- 
ments with weak gravitational lensing offers, in principle, 
a path to unbounded precision in the growth rate / from 
a survey of finite volume in the Universe. In addition such 
techniques allow determination of / without any degener- 
acy with galaxy bias or the power spectrum, i.e. one can 
measure /, not just f /b or f^P. When the redshift sur- 
vey makes use of dark-matter halos (or the galaxies within 



them) , a practical survey sees very shallow gains in precision 
on / with increased redshift-survey target density, because 
only relatively rare halos have the substantial range in bias 
that is necessary to realize large gains from the MS bias- 
modulation technique. Nonetheless we show that the addi- 
tion of lensing information to redshift surveys yields very 
large improvements in precision on growth parameters of 
all kinds. Typically the precision available from the com- 
bined survey technique is equivalent to a 10-fold increase in 
volume of a redshift-only survey. In other words, the con- 
straints available from the combined survey would require 
10 observable Universes without the lensing information on 
biases. 

This significant improvement increases the potential for 
growth tests at low {z < 1) redshifts since concerns about 
the small observable volume are less important. Low-redshift 
galaxy surveys have several advantages, the most obvious 
being higher fluxes hence less expensive redshift surveys, 
particularly if we have techniques that benefit from nP ^ 1. 
Another advantage to z < 1 redshift surveys is that one can 
image and resolve a large number of galaxies that are lensed 
by the spectroscopic targets, i.e. are in the background. At 
z < 1, massive halos are helpfully marked by well-developed 
red sequences. Our naive expectation is that the most bla- 
tant violations of GR should be found in the acceleration 
epoch. And the techniques herein ameliorate the greatest 
disadvantage of 2; < 1 surveys, which is the limited observ- 
able volume. 

Several issues must be addressed before the growth tests 
forecast herein could be realized. First, the Kaiser formulae 
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do not hold to percent accuracy in A'^-body tests and im- 
proved predictions for non-linear behavior in redshift space 
are needed — this is true for all the redshift-space distortion 
methods, not just the combined lensing-|-redshift surveys. 
Particularly crucial for the combined method is that the co- 
variance in Equation (|1U|I between galaxy bins be accurate, 
namely that the velocity field is unbiased and the real-space 
covariance between galaxy bins be the product of each bin's 
covariance with the mass. Other issues we will address in 
future work are the impact of projections and shape noise 
on the determination of bias via weak lensing. The optimal 
depth for the redshift survey is also an issue: Figure indi- 
cates limited gains for surveys of halos < 10^^h~^MQ , but it 
is possible that the determination of bias by lensing-galaxy 
cross-correlation will require a deeper galaxy sample in order 
to reach low levels of stochasticity and exploit the absence 
of sample variance. 

This work is supported by grant AST-0908027 from the 
NSF and DOE grant DE-FG02-95ER40893. 
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